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1 Introduction 



During a century the only cases of integrability of the Euler-Poisson equations were the isotropic case 
and the cases of Euler (1758) and Lagrange (1788). In 1888 S. Kowalevski found a new highly non- 
trivial case of integrability [8]. In modern terms, this is an integrable system on the e(3) algebra with 
quadratic and quartic (in angular momenta) integrals of motion. 

Furthermore, by using a mysterious change of variables, she showed that equations of motion 
for the new case of integrability are linearized on the abelian variety by means of the Jacobi-Abel 
theorem about the inversion of a system of abelian integrals [8]. At the moment no separation which 
is alternative to her original separation of variables is known for this system , even though there is a 
large body of literature dedicated to the problem, including the detailed geometric description of the 
invariant surfaces on which the motion evolves, see books [TJ [5] and references within. 

In this paper we discuss the direct method of finding variables of separation without any additional 
information (ingenious and at times obscure change of variables, Lax matrices, r-matrices, links with 
soliton equations etc). For example, we apply the machinery of bi-Hamiltonian geometry to the 
Kowalevski top at zero level of the cyclic integral of motion, which is a particular case of the generic 
Kowalevski top. The rational Poisson bivector associated with famous Kowalevski variables may be 
found in [53]. Here we specially do not consider Kowalevski variables in order to get only the new 
variables of separation and the new underlying polynomial Poisson structures. 

The other aim is the construction of different variables of separation lying on the distinct algebraic 
curves [26] . Relations between such distinct curves give us a lot of new examples of reductions of 
Abelian integrals and, therefore, they may be the source of new ideas in the number theory, algebraic 
geometry and modern cryptography [TJ |3J [TU] . 

In Section 2 we construct new compatible Poisson bivectors for the Kowalevski top. In Section 
3 we find the new corresponding variables of separation and the separated relations. Finally, some 
concluding remarks can be found in the last Section. 
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2 The bi-hamiltonian structure 



2.1 Description of the model 

According to 8] , the Kowalevski top is a dynamical system with the following integrals of motion 



Hi 



Ji + J. 



2Jl 



CiXi, 



ci e 



(2.1) 



H 2 = {J( + J^Y -2[x 1 {Ji- J z 2 ) + 2x 2 JiJ2)ci + {x{ +xi)c{. 

Here Ji are the components of the angular momentum in the moving frame of coordinates attached to 
the principal axes of inertia. The position of a rigid body is fixed by the components Xi of the Poisson 
vector, which are the cosines between the axes of the body frame and the field up to a constant. 

Using the Hamilton function Hi and the Lie-Poisson bracket {., .} on the Euclidean algebra e*(3) 
the customary Euler-Poisson equations may be rewritten in the hamiltonian form 



Ji — { Ji, Hi} , 



= {x u Hi} , where {/, g} = (Pdf, dg) 



(2.2) 



In coordinates z — (xi, X2, £3, Ji, J2, J3) on e*(3) the Lie-Poisson bivector P is the following antisym- 
metric matrix 
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It has two Casimir elements 

Pda 2 = 0, 



Ci 



fc=i 



Co 



(x,j) = y^x k Jk- 



(2.3) 



k=l 



After fixing values of the Casimir elements 

Ci 



one gets a generic four-dimensional symplectic leaf O a b, which is topologically equivalent to the cotan- 
gent bundle T*S 2 of the sphere S 2 with radius a. However, the symplectic structure of O a b is different 
from the standard symplectic structure on T*S 2 by the magnetic term proportional to b [13] . 

The Kowalevski top is an integrable system on the phase space O a b because the two independent 
integrals of motion -Hi. 2 (|2.1|) are in the involution 



{Hi,H 2 }=(PdH 1 ,dH 2 )=0. 



(2.4) 



In mechanics the Casimir function Ci is a norm of the unit Poisson vector such as a — 1, whereas 
second Casimir function C2 is called a square or cyclic integral of motion [2j [27] . 

Remark 1 In original Kowalevski work the first step in the separation of variables method is the 
complcxification: she considers 

Zl = Ji + iJ2, z 2 = J\ — 1J2 
as independent complex variables. Next she makes her famous change of variables 



Sl,2 



fl(zi,Z 2 ) ± >/fl(zi,Zi)fi(z2,Z2) 



2( Zl -z 2 ) 2 

The fourth degree polynomials R(zi, z^) we will not specify here. It brings the system 

(-l) fe (si -s 2 )s k = ^P(s k ) , 



to the form 



1,2, 



2 



where 



P(s) =A[(s- H 



K 



2„2 



s((s- H) A + c{a 



K 



c\b 



(2.5) 



Here the pairs ^Sfc, T)k — \J P(sk)j , can be regarded as coordinates of points on the Kowalevski curve 
of genus two 

C kow : V 2 - P(s) = . (2.6) 
At zero level of the cyclic integral of motion C2 = the Kowalevski curve has the same genus 2. 

We address the problem of separation of variables for the Hamilton- Jacobi equation as well. At 
C2 = the symplectic leaves O a o are completely symplectomorphic to T*S 2 [13]. We will only consider 
such symplectic leaves and, therefore, all the formulae below hold true up to C2 = 0. 

For this Kowalevski top on the two-dimensional sphere we want to calculate different variables of 
separation and, according to the general usage of the bi-hamiltonian geometry, firstly we have to find 
the second dynamical Poisson bivector P' equipped with some necessary properties [5] [9] . 



2.2 Dynamical Poisson bivectors 

According to [HI [22J, [23l [25] let us suppose that the desired second Poisson bivector P' is the Lie 
derivative of P along some unknown Liouville vector field X 

P' = C X {P). (2.7) 

In addition it has to satisfy the following equations 

[P>,P'] = [£ x (P),£x(P)]=0, (2.8) 

and 

{H lt Ha}' = (P'dH u dH 2 ) = 0, (2.9) 

where [., .] is the Schouten bracket. 

The first assumption (|2.7p guarantees that this dynamical bivector P' is compatible with the 
given kinematic Poisson bivector P, i.e. [P, P'\ = 0. In geometry such bivector P' is said to be the 
2-coboundary associated with the Liouville vector field X in the Poisson-Lichnerowicz cohomology 
defined by P. 

The second condition l[2.8jl means that P' is the Poisson bivector, i.e. that the Jacobi identity 
is true. The third equation (|2.9[) relates P' with the given integrable system. In the wake of this 
agreement the foliation defined by the H12 is the bi-Lagrangian foliation [5[[9]. 

The system of equations (j2.8H2.9p has infinitely many solutions with respect to X (20] [23] . So, in 
order to get some particular solution we have to narrow the search space. In this paper we suppose 
that 

P'dd, 2 = 0, (2.10) 

and that the components Xj of the Liouville vector field X — ^2 Xj dj are non-homogeneous polyno- 
mials in momenta Jk 

N m 
m=0 k=0 

with unknown coefficients g(xi, ^2^3) [21, 22, 25]. Here we explicitly use the restriction C2 = 0, i.e. 
that J 3 = — (xi J\ + x 2 J2)/x 3 . 

Upon substituting this polynomial ansatze into the equations (|2.8l2.9ti2.10[) and demanding that 
all the coefficients at powers of Jk vanish one gets the over determined system of algebro-differential 
equations. Such systems are solved on personal computer by using modern software in a few seconds. 
So, the only real problem is the classification and the analysis of the received computer results. 
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The first three nontrivial solutions arise only in the cubic case N = 3. Components of the first 
real vector field X^ are equal to 

y/x\ + x\(x 1 J 1 -x 2 J 2 ) J3 (1) V x i+ X 2 (%iJl ~ ^2^2) (1) 

1 2x1X3 J 2x2X3 3 

^ - -x/^fHaP-^+T^^^ (2-11) 



(xf + x 2 ) J 3 


Jf J3 




2X1X3 


(xf + x 2 ) Jf 


J 2 J 3 


2 "^3 


2x 2 x 3 



X« = -y^+x 2 
The components of the second real vector field X^ read as 

(2) 2{x\ + X 2 .) (2 ) 2x X (x 2 + X 2 ) (2 ) _ 

A, — ^1^3; A 2 — J lV3 ) A 3 — 

X 3 " X 2 X 3 

(2 ) _ (xf +X 2 ) 2 3 / X1X3 \ 2 3^+z| Cl X 2 J 2 

* 4 " 3x 2 x§ Jl ~( v Jl + ^xf J3 J J3+ ^T 1 2 3 + ^~ ' 



(2) (x^+x 2 ) 2 t3 / (2x 2 -x 2 )x 3 2(x|+x|) \ 2 

* 5 - 3x 2 x 2 J i-{ J 3- 3x 2 x 2 J3 " xix 2 Jl)Js 
xl + xl j j j ci(2xi J 2 - x 2 Ji) 



J1J2J3 , 

X1X3 z 

(2) _ 2x 2 + x 2 . 3 cix 2 (xi J 2 - X 2 Jl) 
6 " 3x 2 3+ 2x3 

The components of the third vector field X^ are the complex functions on initial variables 

„(3) ix 2 (xi +ix 2 ) 2 2 2x 2 (xi +ix 2 ) • r~T 

Al — o 2 JlJ2, 1— \ — J- , 

Xi Xi 

x ( 3) = i(xi+ix 2 ) 2 j2 _ 2(xi+ia;2)JiJ2; x (3) =0j (2.13) 

Xi 



v (3) , r ; t \ t2 1 t3 , 2(2xi + ^2)^3 3 4x 2 j3 , lX 2 (x 1 J 1 - X 3 J 3 ) 

4 = ^ ~ ^ 2 ~ 3 1 3x| 3 3zi a? J 

+ C1X3J3, 

y(3) 2i 7 3 / j ■ j w j 1 r3 2ia; 3 .3 , ^|(2x! -ix 2 ) 3 ixjxs 2 , 

5 ~ 3" 1 ~ ( ~ ' ~ 3 2 ~ 3^~ 3xf 2 ~ "xl 3 " 2 

+ icix 3 J 3 , 

= \ <3 + ^ - ^ J| - ( J 2 + jf ) J 3 - C^X! + ix 2 ) J 3 ■ 
O Xi 

The quartic ansatze yields a lot of solutions, which will be classified and studied in future. 
Let us show the simplest part of these real and complex Poisson brackets explicitly 



r 1 _ r 1(1) _ \/ x \ ~b x 2 J3 

\Xj, Xjf — &ijk Xk j \ x ii x jJ — &ijk Xk 



X3 

! 1(2) _ ^( x l "b x 2) r -1(3) 



2ieijk (ix 3 J 3 - xi J 2 + x 2 Ji)x k ■ 

Here Sijk is the totally skew-symmetric tensor. Other brackets are appreciably more tedious expres- 
sions. The complex Poisson structure may be rewritten in the lucid form by using the 2x2 Lax 
matrices [7J [T5] and the bi-hamiltonian structure associated with the reflection equation algebra [51] . 

It is easy to prove that the corresponding Poisson bivectors P^\ P^ 2 ' and P^ have the following 
properties 

[P«p( 2 )] =0, [P* 1 ',?^]^, [P( 2 \P^] ^0 (2.14) 
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with respect to the Schouten brackets. It means that P^ and P^ are compatible bivectors, whereas 
the complex bivector P^ is incompatible with the real bivectors. 

Remark 2 For any bivectors P and P' there are a lot of vector fields X, such as P' = Cx(P)- Above 
we put X% = in order to restrict this freedom. It may be the origin of some non-symmetry and 
irregularity in expressions (|2.11I2.12[) and (|2.13[) . 

Remark 3 There are two rational Poisson bivectors P' for the Kowalevski top. The first bivector is 
associated with the Kowalevski variables of separation and the underlying reflection equation algebra 
[24j . The second bivector is related with the Lax matrix of Reyman-Semenov-Tian-Shansky and the 
linear r- matrix algebra [19]. The components of the corresponding vector fields X are logarithmic 
functions in momenta. 

To sum up, using the applicable polynomial ansatze for the Liouville vector field X we got two 
compatible real cubic bivectors P^ 1 ' 2 ' = £^(1,2) (P) and one complex cubic bivector P^ = C X (3) (P) for 
the Kowalevski top. Although these bivectors are defined by arbitrary value of C2, they are compatible 
with the initial Poisson bivector P only for C2 = 0. The application of this Poisson bivectors will be 
given in the next section. 



3 Calculation of the variables of separation and the separation 
relations 

A system of canonical variables (qi, . . . , q n ,Pi, ■ ■ ■ ,p n ) 

{qi,Qk} = {pi,Pk} = 0, {qi,pk} = kk (3.1) 

will be called separated if there are n relations of the form 



3>i(qi,Pi, Hi,...,H n ) =0 , i = l,...,n, with det 



dHi 



t^O, (3.2) 



binding together each pair (qi,Pi) and Hi,..., H n . 

The reason for this definition is that the stationary Hamilton- Jacobi equations for the Hamiltonians 
Hi = cti can be collectively solved by the additively separated complete integral 

n 

W(qi,...,q n ;ai,...,a n ) = }^Wi(q i ;ai,...,a n ) , (3.3) 

i=l 

where Wi are found by the quadratures as the solutions of the ordinary differential equations. 

The integrals of motion (Hi H n ) are the Stackel separable integrals if the separation relations 
(|3.2j) are given by affine equations in Hj, that is, 

n 

/,Sij(qi,Pi)Hj - Uj(qi,pi) = , i = l,...,n, (3.4) 

with an invertible matrix S. The functions Sy and Ui depend only on one pair (qi,Pi) of the canonical 
variables of separation, it means that 

{Sik,qj} = {Sik,pj} = {Sik,Sj m } = 0, i^j, (3.5) 

and similar to U 

{Ui, qj } = {Ui, Pj } = {U u Ui) = 0, i ? j. (3.6) 
In this case S is called the Stackel matrix, and U the Stackel potential. 

Remark 4 We have to point out that the definition of the Stackel separability depends on the choice 
of Hi. Indeed, if (Hi, . . . , H n ) are Stackel-separable, then Hi = Hi(Hi, . . . , H n ) will not, in general, 
fulfill the affine relations of the form p.4[) . 
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Remark 5 The method of the separation of variables for a long time served an important, but techni- 
cal role in solving Liouville integrable systems of classical mechanics. A new, and much more exciting 
application of the method came with the development of quantum integrable systems. Because of 
the fact that the quantization of the action variables seems to be a rather formidable task, quantum 
separation of variables became an inevitable refuge. In fact, it can be successfully performed for many 
families of integrable systems with affine separated relations (|3.4p only [THl [T7] . 

So, our second step is the calculation of the canonical variables of separation (qi,pi) and of the 
separation relations (13. 2j) . According to [3|9], the coordinates of separation m are eigenvalues of the 
recursion operator^ which are the so called Darboux-Nijenhuis variables. In order to get the recursion 
operator TV = P'P -1 we have to find restrictions P, P' of the Poisson bivectors P and P' onto the 
symplectic leaves. 

We can avoid the procedure of restriction using the n x n control matrix F defined by 

n 

P'dH = P(PdH), or F 'dH. t = P^Py dHj, i = l,...,n. (3.7) 

i=i 

The bi-involutivity of the integrals of motion (|2.4I2.9[) is equivalent to the existence of F, whereas 
the imposed condition (|2.10p ensures that F is a non-degenerate matrix. In this case eigenvalues 
of this matrix coincide with the Darboux-Nijenhuis variables and we can easily calculate the desired 
coordinates of separation qi. 

Moreover, for the Stackel separable systems the suitable normalized left eigenvectors of the control 
matrix F form the Stackel matrix S 

F = S~ 1 di8,g{q 1 ,...,q n )S 



which would allow us to get separated relations 

So, the main problems are the finding of the conjugated momenta pi and the construction of the 
separation relations cf>j (|3.2[) for the generic non-Stackel separable systems. Below we show how we 
can solve these problems using the same control matrix F and some additional useful observations. 

3.1 The real compatible Poisson bivectors 

For the first Poisson bivector PW (|2.11[) the entries of the control matrix PW read as 
(i) = (2x 2 + 2x 2 +x 2 )(J 2 + J 2 ) (1) _ 1 

11 . o / 9 o" 12 



^x\\Jx\ + x\ " %y/x\ + X 



2 



r 2l 



(2x1 + 2xj + xl)(Jl + J$) cis/x'f + xi{x x (Jl - Jf ) + 2x 2 Ji J 2 ) c\yjx\+xl 



F 22 



(i) _ J i + J ? 



2x^ y/ X 2 -\- X^ 
2 



2y/x\ + x\ 

The eigenvalues of this matrix are the required variables of separation q\ t 2 
det(F^~\I) = (X- qi )(X-q 2 ) 



2 ^x 2 +x 2 (J 2 + J 2 ) c 1 [2x 1 {J l - J 2 ) + 4X2^1 J 2 + c x x, 

A — — n A — 



The matrix of normalized eigenvectors of P^ 1 ) does not form the Stackel matrix, because property 
(|3.5[) is missed, and the underlying separation relations differ from the Stackel affine equations 
in H\ 2- 
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For the second Poisson bivector P^ 2 ' (12.12j) the entries of the control matrix are equal to 
J? + Jf - cjxi (xi J 2 - x 2 Ji) 2 



F ( 2 ) 



P (2) 



(2) 



= -Vi + JiY 1 + 



2(a;? + z|) 



F, 



+ ci(x 2 + x 2 ) 



1 

4' 



'2(xi{Jl- Ji)+2x 2 JiJi 



Cl 



F 



(2) 



J 2 + Jf + 2 Jf + cixi 



The eigenvalues /i l2 of the matrix F^ 2 ) are the roots of the equation 



det(p( 2 >-A7) = (A-A)(A-/ 2 ) 



= A 2 - c\X\ - 



(x 2 Ji - xiJ 2 - x 3 J 3 )(x 2 J\ - x x J 2 + x 3 J 3 ) 



(2(x 2 J\ - x 1 J 2 )J 3 - c x x 2 x 3 ) 2 
4xj 

Remark 6 According to [5]IH] the compatibility of pf 1 ' 2 ) (|2.14[) ensures that in the Darboux-Nijenhuis 
variables q,p the corresponding restrictions of p( 1,2 ) look like 



p(i) = 



/ qi \ 

q 2 

-qi 

V o -q 2 o o y 



p(2) = 



/ 








h 





\ 













h 




-fx 













V 





-7*2 











where /i l2 are the functions on q,p such as 

tiijj} = {Phfj} = 0) Mi' 
So, /i is the function only on gi and p\ and similar / 2 is the function on q 2 and p 2 . 

We can find these functions j\_ 2 using the Poisson bracket. Namely, it is easy to see that the 
recurrence chain 



<t>X = {fx (Qx,Px), <li}, fa = {<t>x,qx}, 



4>i = 



(3.8) 



breaks down on the third step <f> 3 — 0. It means that f± is the second order polynomial in momenta 
Pi and, therefore, we can define this unknown momenta in the following way 



Pi 



2x 3 ^4(lZI 2 Ji — Xi J 2 ) qi + c\yjx\+ x\J 2 



^ (4 v /x 2 +x 2 (J 1 2 + J 2 ) gi + ci .ti(J 2 - J 2 ) + 2x 2 Ji J- 



(3.9) 



up to canonical transformations p\ — > p\ + g(gi). 

The similar calculation for the function f 2 {q 2 ,p 2 ) yields the definition of the second momenta 



2x 3 [4(x 2 Ji - Xi J 2 ) q 2 + ciyjx\+ x\ J 2 



92 



{4^/x 2 +x 2 {J 2 + Jf ) q 2 + ci (xi( J\ - Jf ) + 2x 2 Ji J 2 ) 



(3.10) 



In fact we have to substitute q 2 instead of q\ only. 

So, one gets the canonical transformation from the initial physical variables (x, J) to the variables 

of separation (q,p) using a pair of compatible bivectors PC 1 ' 2 ) and the corresponding control matrices 
p{x,2)_ 

In these separated variables entries of the matrix S of normalized eigenvectors of F^ depend on 
the pair of variables {qi,Pi) and the Hamilton function 



5,i = -2H X - [4q 2 - 4 p 2 , S 1)2 = 1, 



i = 1,2. 
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In this case S may be called the generalized Stackel matrix. The separated relations (|3.4|) look like 
SnHi +H 2 - [Hi - 14 - - a 2 (c 2 - 16^)) =0, i = 1, 2. 

and, therefore, the generalized Stackel potential C/j depends on the Hamilton function too. 

So, we can say that the variables of separation (qi,Pi) lie on the algebraic hyperelliptic curve C of 
genus three defined by 



- a 2 (c 2 - I6q 2 ) = 0. (3.11) 

This equation is invariant with respect to involution (q,p) — > (—q,p)- Factorization with respect to 
this involution give rise to elliptic curve 



- a 2 {c\ - 16z) = 0, z^g 2 (3.12) 

Due to the standard formalism we have to calculate differential on this curve 

dz 



z(z,py 

Then it's easy to prove that 



Z(z, p) = p [c\ - 16z) Uh-l - p 2 {c\ - 16z)) . 



Qi , 92 _ n (c 2 - 16g 2 )p 2 gi (c 2 - lGq 2 )p 2 2 q 2 1 , 



Z(q 2 , Pl ) Z{qlp 2 ) ' ^( 9 2 ,Pi) 2-(g 2 2 ,P2) 4 

and 

Pi , 



Z{q 2 ,p) J Z(q 2 ,p) 

qi {c\ - 16q 2 )p 2 dq I" 12 {c\ - 16q 2 )p 2 dq _ t 

ZWJ) + J Z{q^p) ~~l + P 2 

where p has to be obtained from (13.111) . 

Remark 7 The equations of motion are linearized on an abelian variety, which is roughly spiking the 
complexified of the corresponding Liouville real torus. So, even though q\^ 2 are the real variables of 
separation we have to solve the Jacobi inversion problem over the complex field, see more detailed 
discussion in [UH]. 

Remark 8 We have to point out the Kowalevski separation of variables leading to hyperelliptic 
quadratures, whereas in the new variables of separation q\ 2 equations of motion are integrable by 
quadratures in terms of elliptic functions. 

The third part of the Jacobi method consists of the construction of new integrable systems starting 
with known variables of separation and some other separated relations [6 . If we substitute our variables 
of separation (q,p) into the following deformation of ()3.1ip 

<S> {d) {p,q) = <5>{p,q)-8d iq -16d 2 q 2 = 0, d u d 2 € R, (3.14) 

we get the following generalization of the initial Hamilton function 

H [ d) = J 2 + J 2 + 2J 3 2 + Cl xl + + % (J 2 + J 2 ) . (3.15) 

Here the main problem is how to get the Hamiltonian to be interesting to physics. For example, in our 
case we obtained the natural Hamiltonian at d 2 = only. For this system only the integrals of motion 
have been known 1271. 
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The apparent problem is that generalized equations p,14|) have not involution (q,p) — > (—q,p) at 
di 0. Thereby equations of motion are related with the hyperelliptic curve of genus three [TTJ [T2] 
instead of elliptic curve. Nevertheless, at d\ ^ 0, d 2 = we have the same equations (|3.13p as above 



<7i 92 



0, 



(3.16) 



Z(qj,Pi) Z(qlp 2 ) 

(^lggjjgfo (cf- lggfjpfgg = _1 
Z(g 2 , Pl ) Z(^,pa) 4 

where pi. 2 satisfy to the deformed equations (|3.14l) . 
3.2 The complex Poisson bivector 

For the cubic in momenta Poisson bivector (|2.13l) the control matrix is equal to 

F (3) = ( 2(jf + Jf + 2J|) + Ci (:ei + ix 2 ) ™ \ 
V 2(J? + J 2 2 ) 2 -2ci(a;i-ix 2 )(Ji+iJ 2 ) 2 / 
and the Darboux-Nijenhuis coordinates A1.2 are the roots of the characteristic polynomial 

p(3) 

det(F - A/) = (A - Ai)(A - A 2 ) = A 2 - F£f\ + -|- . (3.17) 

As above we can get the conjugated momenta /ix,2 by using compatible with P^ 3 ) bivector of fourth 
order in momenta However we can do it without such calculations as well. 

It is easy to see that in this case matrix S of normalized eigenvectors of is the standard 
Stackel matrix 

Sii = — 2Aj, Si t 2 = lj £=1,2, 
and, therefore, the Stackel potentials 

Ui, 2 = -(S il H 1 + H 2 ) 

are some functions on (Ax,/tti) and (A2,/X2), respectively. 

In fact notion of the Stackel potentials allows us to find the unknown conjugated momenta fii_ 2 
using the Poisson brackets only. Namely, the following recurrence chain of the Poisson brackets 

<k ={Ai,t/x}, <f> 2 = {Ai,<M,..., ^{Ai,^-!} (3.18) 

is a quasi-periodic chain 

03 = 16Ai0i. 

It means that the Stackel potential 17% is a trigonometric function on momenta [i\ and, therefore, we 
can determine this desired momenta 



= <ys(Ai)ln(Vl6Ai0i + 2 ) 



up to canonical transformations ^x — > Ml + sCVO- Here the function <p(\i) is easily calculated from 
{Ax./ux} = 1. 

These variables of separation (Ai,/ii) lie on the hyperelliptic curve of genus three 

C: $(A,^) =e 4iV ^ + ^|L- 4iv ^ I + A 2 - 2H 1 X + H 2 =0. (3.19) 

According to [10j . this curve C are related with an elliptic curve £ and equations of motion for the 
Kowalevski top are linearized on the corresponding abelian variety. 

One main difference is that the variables of separation Ai i2 are complex functions on the initial 
variables {x, J), whereas q\ i2 are real functions on them. It will be important when we express initial 
real variables via real or complex variables of separation after solving of the Jacobi inversion problem 
over the complex field. 

The other difference is that the affine relations of separations (|3.19[) allows us to study quantum 
counterpart of the Kowalevski top [17] . For real variables of separation the procedure of quantization 
is unknown. 
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Remark 9 In framework of the Sklyanin formalism [16] variables of separation are the poles of the 
Baker- Akhiezer function with suitable normalization. In [TS] we find such variables of separation 
mi, 2 for the Kowalevski-Goryachev-Chaplygin gyrostat 

Hi = J\ + Jf + 2J 3 2 + p,J 3 + cixi + c 2 (x\ - x\) + c 3 x lX2 + % (3.20) 

x 3 

using 2x2 Lax matrix, its Baker- Akhiezer vector-function and the reflection equation algebra. 

It is easy to prove that the Darboux-Nijenhuis variables Ai^ Q3-17P are related with the poles u 12 
of the Baker- Akhiezer function by the following point transformation 

Ai,2 = < 2 , (3.21) 
which gives rise to a ramified two-sheeted covering of C, see [TO] , 

4 Conclusion 

Starting with the integrals of motion for the Kowalevski top we found three polynomial in momenta 
Poisson bivectors, which are compatible with the canonical Poisson bivector on the cotangent bundle 
T*S 2 of two-dimensional sphere. 

Then in framework of the bi-hamiltonian geometry we get new real variables of separation (q,p) 
for the Kowalevski top on the sphere and reproduce known complex variables (X,fi). These variables 
are related by the canonical transformation 

Al,2 = Al,2(<?l,<72,Pl,P2), Ml.2 = Ml,2(<Zl,?2,.Pl,P2), 

which may be rewritten as a quasi-point canonical transformation [14] 

Ai,2 = Ai ! 2(9ij 92, Hi,H 2 ), 

which relates two hyperelliptic curves of genus three. We can assume that it is no rational cover and 
that these Jacobians are non-isogeneous in Richelot sense [15]. Similar transformations relate these 
curves with the Kowalevski curve of genus two. Further inquiry of such relations between hyperelliptic 
curves goes beyond the scope of this paper, see discussion in [TJ 131 [TU]. 

The proposed approach may be useful for the investigation of other integrable systems with inte- 
grals of motion higher order in momenta, for instance, the search of another real variables of separation 
for the Kowalevski-Goryachev-Chaplygin gyrostat and its various generalizations [27] . 
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